library(policyPlot)

get_fofnet <- function(message = TRUE) { 
  #read path from text file
  path_file <- "[insert path to public FoF data]"
  
}

dfa_splits <- function(group = "networth", agg_level = "default", subgroup = "", condition_cat = NULL) {
  if (is.null(condition_cat)) {
    if (group == "networth") {
      if (agg_level == "bus_cycle2") {
        groupnames <- c("TopPt1", "RemainingTop1", "pct90to99", "pct70to90", "pct50to70", "Bottom50")
        scf.recode <- c(`1` = "Bottom50", `2` = "Bottom50", `3` = "Bottom50", `4` = "Bottom50", `5` = "Bottom50",
                        `6` = "pct50to70", `7` = "pct50to70", `8` = "pct70to90", `9` = "pct70to90", `10` = "pct90to99",
                        `11` = "pct90to99", `12` = "RemainingTop1", `13` = "TopPt1")
        tag <- "scf_fa_nw"
        input.file <- "dfa_scf_fa_nw_buscycle_detail"
        agg_vec <- c(.1,.9, 9, 20, 20, 50)
      }
    }
  }
  if (exists("agg_vec") == FALSE) {
    output <- list(groupnames, scf.recode, input.file, tag)
    names(output) <- list("groupnames", "scf.recode", "input.file", "tag")
  }
  
  if (exists("agg_vec") == TRUE) {
    output <- list(groupnames, scf.recode, input.file, tag, agg_vec)
    names(output) <- list("groupnames", "scf.recode", "input.file", "tag", "agg_vec")
  }
  
  
  return(output)
}


## Write a function that brings in SCF output
scf_balance_sheet <- function(excel, hh = FALSE, income = FALSE){
  
  sheet_names <- excel_sheets(excel)
  sheet_names <- sapply(strsplit(sheet_names," "), `[`, 1)
  
  ## Read in the excel files. The format must be consistent (check DFA_reconciliation). 4 sheets: Nonfinancial Assets, Financial Assets, Other Assets, and Liabilites. 
  ###!!! Need more extensive documentation on how the SCF inputs need to be formated 
  non_fin_sheet <- read_excel(excel, sheet = which(sheet_names == "nonfin"), trim_ws = TRUE, col_types = "numeric") %>%
    drop_na(2)
  fin_sheet <- read_excel(excel, sheet = which(sheet_names == "fin"), trim_ws = TRUE, col_types = "numeric") %>%
    drop_na(2)
  other_sheet <- read_excel(excel, sheet = which(sheet_names == "othassets"), trim_ws = TRUE, col_types = "numeric") %>%
    drop_na(2)
  liabilities_sheet <- read_excel(excel, sheet = which(sheet_names == "liab"), trim_ws = TRUE, col_types = "numeric") %>%
    drop_na(2)
  
  if (hh == TRUE) {
    aux_sheet <- read_excel(excel, sheet = which(sheet_names == "households"), trim_ws = TRUE, col_types = c("numeric", "numeric", "numeric", "skip", "skip", "skip", "skip"), .name_repair = "minimal") %>%
      drop_na(2)
  }
  
  balance_sheet <- list()
  
  balance_sheet[["real_estate"]] <- non_fin_sheet$scf_fa_real_estate
  balance_sheet[["con_durables"]] <- non_fin_sheet$scf_fa_consumer_durables
  
  balance_sheet[["check_curr"]] <- fin_sheet$scf_fa_check_fgn_dep_curr
  balance_sheet[["time_dep"]] <- fin_sheet$scf_fa_time_dep
  balance_sheet[["mmfs"]] <- fin_sheet$scf_fa_mmmf_shares
  
  balance_sheet[["us_gov_muni"]] <- fin_sheet$scf_fa_gov_muni_bnds
  balance_sheet[["corp_for_bonds"]] <- fin_sheet$scf_fa_corp_for_bnds
  
  balance_sheet[["oth_loans"]] <- other_sheet$scf_fa_othln_owed
  balance_sheet[["morts_owed"]] <- other_sheet$scf_fa_mort_owed
  
  balance_sheet[["corp_equ_mfs"]] <- other_sheet$scf_fa_corpeq_mut
  balance_sheet[["life_ins"]] <- other_sheet$scf_fa_life_ins
  balance_sheet[["pension"]] <- other_sheet$scf_fa_pension
  balance_sheet[["non_corp_eq"]] <- other_sheet$scf_fa_equity_non_corp
  balance_sheet[["misc_assets"]] <- other_sheet$scf_fa_misc_assets
  
  balance_sheet[["morts"]] <- liabilities_sheet$scf_fa_mort_debt
  balance_sheet[["cons_cred"]] <- liabilities_sheet$scf_fa_consumer_credit
  balance_sheet[["dep_loans"]] <- liabilities_sheet$scf_fa_dep_inst_loans
  balance_sheet[["oth_loans_adv"]] <- liabilities_sheet$scf_fa_oth_loans
  
  balance_sheet[["unpaid_life_ins"]] <- liabilities_sheet$scf_fa_life_ins_unpaid
  
  if (hh == TRUE) {
    balance_sheet[["hhcount"]] <- aux_sheet$hhcount
  }
  
  if (income == TRUE) {
    
    # use whichever income data are available -- varies by split
    
    if ("inc" %in% sheet_names) {
      
      inc_sheet <- read_excel(excel, paste0(sheet = "inc by ", tag), col_types = "numeric") %>%
        drop_na()
      balance_sheet[["inc"]] <- inc_sheet[["income"]]
    }
    
    if ("nipainc" %in% sheet_names) {
      nipainc_sheet <- read_excel(excel, sheet = paste0("nipainc by ", tag), col_types = "numeric") %>%
        drop_na()
      balance_sheet[["nipainc"]] <- nipainc_sheet[["nipa_equiv_income"]]
    }
    
    if ("norminc" %in% sheet_names) {
      norminc_sheet <- read_excel(excel, sheet = paste0("norminc by ", tag), col_types = "numeric") %>%
        drop_na()
      balance_sheet[["norminc"]] <- norminc_sheet[["norminc"]]
    }
  }
  
  balance_sheet <- bind_cols(balance_sheet)
  
  ## Now construct the balance sheet in the format of the B101h table. The series labels are done manually to look as much like the fame tables as possible.
  
  ids <- non_fin_sheet[,1:2] %>%
    mutate(date = as.yearqtr(year + .5)) %>%
    setNames(c("year", "cat", "date")) %>%
    select(date, cat)
  
  
  scf_data <- bind_cols(ids, balance_sheet)
  
  return(scf_data)
}

dfa_getfame <- function(series, database) {
  df <- policyPlot::fame2df(series, database) %>%
    drop_na() %>%
    mutate(date = as.yearqtr(date)) %>%
    mutate_if(is.array, as.numeric)%>%
    dplyr::filter(date >= as.yearqtr("1989q3"), date <= as.yearqtr(lastqtr)) %>% # _filter out observations before the start date (1989q3)
    group_by(date) %>%
    summarise_all(data.table::last) # convert to quarterly series by choosing last value in each quarter
  
  return(df)
}


dfa_file_names <- function(item = "", extension = ".csv", group = "networth", method = "fernandez", agg_level = "default", scf_ignore = NULL, filetag = NULL) {
  
  name_vec <- character()
  name_vec[1] <- "dfa"
  name_vec[2] <- group
  index <- 3
  
  if (method != "fernandez") {
    name_vec[index] <- method
    index <- index + 1
  }
  
  # if (agg_level != "default") {
  #   name_vec[index] <- agg_level
  #   index <- index + 1 
  # }
  
  if (!is.null(scf_ignore)) {
    name_vec[index] <- paste0("skip_", scf_ignore) 
    index <- index + 1
  }
  
  if (item != "") {
    name_vec[index] <- item
    index <- index + 1
  }
  
  if (!is.null(filetag)) {
    name_vec[index] <- filetag
  }
  
  file_name <- paste(name_vec, collapse = "_")
  file_name <- paste0(file_name, extension)
  
  return(file_name)
}

scale2one <- function(vec) {return(vec/max(na.omit(vec)))}

neg_to_zero <- function(x) ifelse(x < 0, 0, x)
